Phase transition of one dimensional bosons with strong disorder 
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We study one-dimensional disordered bosons at large commensurate filling. Using a real-space 
renormalization group approach, we find a new random fixed point which controls a phase transition 
from a superfluid to an incompressible Mott-glass. The transition can be tuned by changing the 
disorder distribution even with vanishing interactions. We derive the properties of the transition, 
which suggest that it is in the Kosterlitz-Thouless universality class, and discuss its physical origin. 



Models of interacting bosons subject to quenched dis- 
order are important for understanding many real sys- 
tems and present a considerable theoretical challenge. 
Bosons in a disordered potential, for instance, were stud- 
ied extensively [1-7] in relation to granular supercon- 
ductors, Josephson-junction arrays, and '^He on porous 
media. Recently, there has been interest in this problem 
in the context of experiments with ultra cold atoms [8, 9]; 
These may provide a controllable environment for study- 
ing quenched disorder and interactions. In this letter we 
take a new perspective, and utilize the real-space renor- 
malization group (RSRG) method [10-12]. Using this 
method it was shown that random spin chains often ex- 
hibit a universal behavior independent of the disorder 
realization. Using RSRG we describe the properties of a 
superfluid-insulator transition of the random 0{2) quan- 
tum rotor model, and find that it exhibits a similar degree 
of universality. 

A standard model of granular superconductors and 
Josephson arrays in 1-d is the 0(2) quantum rotor Hamil- 
tonian: 

1 f d \^ 

H = -J2uJ-i—) -J2jjCos{ipj+i-ipj). (1) 

where Ui is the grain charging energy, and Ji is the 
Josephson coupling between grains. For uniform cou- 
plings, the zero-temperature quantum partition function 
of the chain can be mapped to that of the two dimen- 
sional classical x-y model with unit coupling and effective 
temperature T'-ff = y/U/J [13]. It follows that the pure 
system displays a Superfluid-insulator transition of the 
Kosterlitz-Thouless (KT) universality class at a critical 
interaction strength \JU j J ^ 1. 

The effect of disorder on (1) is important for under- 
standing real systems. In systems such as granular su- 
perconductors, fluctuations in grain size, and distance be- 
tween grains are the main source of randomness, which is 
exhibited in J and U. The analysis of the random model 
is often carried out using the uniform model because ran- 
domness in U and J is perturbatively irrelevant. But 
the irrelevance of weak disorder does not rule out fixed 
points, and hence possible new phase transitions, at finite 
disorder. We add that diagonal (potential) disorder can 



become relevant above a critical interaction strength [1] . 
Here assume that our model is strictly particle-hole sym- 
metric, and omit diagonal disorder form (1). 

In this letter we identify a new quantum phase tran- 
sition from a disordered superfluid to an incompressible 
gapless phase (Mott-glass) using the RSRG method. Re- 
markably, the fixed point that controls the transition cor- 
responds to a classical model {Ui =0). It can thus be 
tuned by varying the disorder in J even at arbitrarily 
small charging energy U, in sharp contrast with the uni- 
form case. We study the universal properties of the tran- 
sition using the RG and show that it is KT like with a 
dynamical exponent z = 1, as in the uniform case. In ad- 
dition to the RG treatment we present a simple physical 
picture of the emerging phases and the phase transition. 
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FIG. 1: Possible fates of the chain after repeated application 
of the RG decimation rules, (a) In the superfiuid phase sites 
are joined into ever growing clusters, (b) Insulating phase, 
where clusters become effectively disconnected at low energy 
scales. 

The Hamiltonian (1) is equivalent to the Bose-Hubbard 
model at large commensurate filling [14], but it differs 
in some respects from the traditional dirty boson prob- 
lem [1, 3, 5]. First, the Hamiltonian (1) obeys a lo- 
cal particle hole symmetry. Thus, it describes only the 
Superfluid-insulator transition at commensurate filling. 
In addition, it does not contain the diagonal disorder 
ivj d/dtpj. The significance of this difference was illu- 
minated by recent Monte-Carlo simulations of large two 
dimensional systems, which considered both types of dis- 
order [15]. While the diagonal disorder is expected to 
ultimately dominate, its effects may be unobservable at 
practical temperatures or length scales. 
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Let us begin our analysis by describing the RSRG 
method as apphed to the random-boson problem. In 
the spirit of Ref. [12], we define the cutoff energy scale 

= maxi {Ui, Ji}. At each RG step, the chain clement 
with largest energy is decimated, renormalizing nearby 
interaction constants. Two types of steps are possible. 
If f7 = Ui for some site i, this site is assumed frozen 
in its lowest energy charge state. A connection between 
nearest neighbors of the site is generated perturbatively: 

Ji-l,i+l = Ji-lJi/^- (2) 

In the alternative case, when f2 = Jj the relative phase 
(fi+i — <fi is assumed frozen in its ground state and the 
two sites are renormalized into a single one. The charging 
energy of the new site corresponds to additive recursion 
relations in terms of capacitances Cj = Cl/Ui: 

Ci=Ci + Ci+i, (3) 



while the hopping connecting the new site with its neigh- 
bors is unchanged. 

Repeated decimations gradually reduce the energy 
scale from its initial value Oq to a lower energy scale 
CI. Depending on the initial distributions of the cou- 
plings, two scenarios emerge. In the first, sites are joined 
together into ever growing superfluid clusters (see Fig. 
1(a)) leading to a superfluid phase. In the second, a grow- 
ing number of sites are eliminated and form disconnected 
clusters leading to an insulating phase (Fig. 1(b)). 

It is convenient to describe the progression of the RG 
transformations and the cutoff energy scale with the pa- 
rameter F = log(rio/^^)- In addition, we define the di- 
mcnsionless coupling Q = d — 1 and /3j = log(f2/Jj), 
characterized by probability distributions f{C,T) and 
g{(3, F) [19]. Their flow with the decreasing energy scale 
is given by the master equations: 



(1 + C)^^^^ + 9o(r) I dCirfC2/(Ci,r)/(C2,r)5(Ci + C2 + 1 - C) + /(C,r) (/o(r) + 1 - go{T)) ., 



dr 



dp 



/o(F) / d;3id/325(/3i,r)5(/32,r)5(/3i+/32-/3)+5(/3,r)(5o(r)-/o(F)), 

I 



(4) 



where goiV) = g{0,r) and /o(r) = /(0,r). The first 

term in each equation describes the flow of the distri- 
bution due to a redefinition of the cutoff following an 
elimination of large couplings. The second term im- 
plements the recursion relations, taking account of the 
renormalized couplings. The last term removes the cou- 
plings neighboring the decimated ones and takes care of 
normalization of the distributions. 

When either typical ^ 1 or /3 ^ 1 one can use 
-^(Ci + C2 + 1 - C) « '^(Ci + C2 - C) in (4). Then the master 
equations are solved by the ansatz: 

/(C, F) = /o(F)e-C/o(r), r) = go{T)e-^^o{T) ^ (5) 
where /o and go obey: 

W = /o(r) - /o(r)ffo(r), ^ = -5o(r)/o(F). 

(6) 

Thus, 

/o(F)=5o(r)-lnffo(r) + A (7) 

It is interesting to note that Eqs. (6) acquire precisely 

the form of the KT flow equations [16] when written in 
terms of y = V/o and x = go. The integration constant 
A controls the flow as depicted in Fig. 2. Distributions 
with A > ~l flow to /o 00 and go ^ 0. This cor- 
responds to vanishing Josephson coupling and U ^ fl, 



namely an insulating phase. For A < —1 and 50 > 1 the 

flow approaches a line of fixed points with a non univer- 
sal 50(00) = 1 -|- a, a > and /o(oo) = 0. This corre- 
sponds to an array of Josephson junctions with no charg- 
ing energy and a power-law distribution of couplings 
p{J) oc (J/Cl)". The critical flow occurs when A = — 1. 
It terminates at the unstable fixed point g* = go{oo) = 1 
(a = 0) and /* = fo{oo) = 0. Note that even a system 
with vanishingly small charging energy and random tun- 
neling can be tuned through the critical point by chang- 
ing the disorder distribution. When a < the system 
flows to the insulating phase. 

We now use Eqs. (6) together with the asymptotic 
forms of go(r), ,/o(r) to show that the transition has KT- 
like universal properties and dynamical exponent z = 
1. First, we establish a connection between energy and 
length scales. As the RG flow proceeds, more sites are 
either eliminated or joined into larger superfluid clusters. 
One can see that the total number of sites at scale F 
decreases as 

Ar(F) = Noe--fo'^f°+3o)dT\ (8) 

The length scale associated with F is the size of super- 
fluid clusters (or bonds between them). In units of the 
original lattice constant it is given by l{T)=No/N{r) = 
exp(/jj"(/o -|- go)dT'). Solving the flow Eqs. (6) near the 
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FIG. 2: Phase space flow of the functions /o(r) and 3o(r). 
Different linos correspond to different values of A in (7). A > 
— 1 (dashed lines) result in flow to the insulating phase. A < 
— 1 (dash-dotted lines) flows either to the insulating phase for 
small go, or to the disordered superfluid phase for large go- 
The solid line ^ = — 1 is the separatrix. The line /o s is a 
line of fixed points which are stable for go > 1 and unstable 
for go < 1. 
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FIG. 3: Numerical RG flows for systems with initial box 
distributions. (a) Distribution of U is centered around 
0.7 + S/2 and that of J around 0.7 - 6/2 with 5 = 
0.03, 0, -0.03, -0.06, -0.09, -0.12 . The width of both dis- 
tributions is always taken to bo 1. /o and go are defined as 
the inverse of the numerical avorago of and /3 respectively. 
The dashed line marks the supcrfluid-iusulator separatrix ac- 
cording to Eq. (7) with ^ = 0. (b) A semilog plot of the dis- 
tributions of the couplings C = ^/U — 1 (top) and = In fJ/ J 
(bottom) at a typical point in the fiow not far from the critical 
point (marked by a black dot in (a)). 



fixed point, = 1, /g = we identify typical time and 
length scales, which characterize the critical flow: 

Tcoce-i/v^, ^cOceVv^, (9) 

where e=l + A>Ois the tuning parameter. The 
exponential divergence toward the critical point suggests 
a transition of the KT universality class with a dynamical 
exponent z = 1. 

Using the universal distributions we found above wc 
can calculate several obscrvables. For example, to find 
the compressibility, we apply the RG to a finite chain 
down to the scale Fi, where only a single site is left 
(iV(ri) — 1). The compressibility extracted from the 
renormalized Hamiltonian of the remaining site Hi = 
U'n? then reads: 

(10) 

Using Eqs. (8), and (6) to derive the asymptotic forms 

of /o(r) and ga{T), wc find in the insulating regime. 
Kin oc (lniVo)/iVo, which vanishes in the thermodynamic 
limit. Similarly, in the SF side Kgf = k/Clo, where k 
is a positive non-imiversal constant, that depends on 
initial distributions. It follows that k is discontinuous 
at the critical point. The gap in the SF phase obeys 
A = fl(Fi)/o(ri)(-ln/o(ri)) ~ y^e'-^. In the insu- 
lating phase, the gap is A ~ n(Fi) ~ 1/lnA'^o- The gap 
vanishes in this case because of rare large clusters that 
have arbitrarily small charging energy. Thus, we identify 
the insulating phase with a Mott glass. 

We have implemented the RG scheme numerically to 
verify that the solution (5) is an attractor of the flow 
for generic initial disorder. We started from chains with 



2.5 X 10^ sites, and different forms of initial disorder dis- 
tributions. The general features of the flow diagram seem 
insensitive to the type of distributions used. A specific 
example calculated from various box distributions of J 
and U is given in Fig. 3. The result clearly shows that 
the distributions converge to the solutions (5), demon- 
strating that they are stable attractors of the RG flow. 
The positions of the minima marked by the dashed line 
in Fig. 3 agree with those predicted by the first correc- 
tion to (6) from expansion of the (^-function in Eq. (4) to 
lowest order in 1/^. 

The validity of the recursion relations (2) and (3) re- 
lies on the distribution of J's being wide in the sense 
that (J) << ri. Clearly this does not hold at the crit- 
ical point where (J) — Vl/2. While the relation (2) can 
be controlled by the smallness of {U) near the critical 
point, corrections to (3) produce next nearest neighbor 
Josephson couplings of order {J)'^ /il. To justify the RG 
scheme used here, we need to show that the longer range 
terms are irrelevant. The corresponding analysis is rather 
complicated and we defer it to a future work. Instead, 
we rederive the phase transition and explain its nature 
using an alternative physical argument. 

To understand the origin of the fixed point distribu- 
tions we employ two simple toy models (see also Ref. 7) . 
First, consider the insulating fixed point. We assume 
that the Josephson couplings take the values Ji = 1 
with probability q, and Ji = with probability 1 — q. 
We also assume that the charging {/, = u <i; 1 is uni- 
form. Such a chain consists of disconnected clusters and 
can describe only the insulating phase. The cluster sizes 
are distributed according to P{N) = , and the effec- 
tive interaction is given by U{N) = u/N. From this 
we immediately see that f{() is indeed given by (5) with 
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/o = I lug I . If we now start eliminating sites with large U 
from the chain, we will exactly reproduce the flow equa- 
tion in (6). Clearly, the insulating phase is stable towards 
weak Josephson couplings between the clusters. There- 
fore the disconnected chain used in this toy model is fairly 
generic once we deal with an insulator. 

We can describe the superfluid fixed line by a chain 
with uniform Josephson couplings Ji = e << 1 and 
charging energies with a binary distribution: Ut ~ 1 
with probability q, and Ui = with probability 1 — q. 
We use perturbation theory to eliminate the sites with 
Ui — 1. Two remaining sites, separated by a block of 
N eliminated sites are effectively connected by J = 
(see Fig. 4). In addition, the size distribution of blocks 
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FIG. 4: Schematic representation of the second toy model. 
The circles on the top graph denote sites with a large charg- 
ing energy (7 = 1. After eliminating these sites we arrive to a 
classical chain with disordered couplings between the neigh- 
boring sites (bottom plot). 

with f/i = 1 is P{N) ~ q^. The syste m is thus described 
by an effective classical chain of Josephson junctions (no 
charging energy) with a power-law coupling distribution: 

PiJ) = P{N)^ = ia + l)iJ/er, (11) 

where a — Inq/ Ine — 1. 

The toy model suggests that the classical chain with 
a > is also a stable fixed line due to the vanish- 
ing probability of finding a J = on the chain. A 
small U perturbation, however, destabilizes the line for 
a < 0. Physically, this is due to vanishing thermody- 
namic stiffness; we prove this in the following. The stiff- 
ness p{L) = L{E{0) — E{0))/9'^ of a given realization of 
a chain of length L is given by p{L) = L/{J2i ^/ Ji)- In 
the limit L — > oo, the average over the chain appear- 
ing in the denominator can be replaced by an average 
over the distribution yielding p(oo) = max{Q, a/(a-|-l)}. 
Clearly, a state with zero stiffness is expected to be un- 
stable to adding a small charging energy Aq. To see 
this more explicitly, consider perturbing the system with 
a small uniform Ui = Aq. Wc choose an energy scale 
Aq < <^ flo — e. The weak bonds with J < il are 
separated by clusters of average size L — {Qq/Q)'^^^ con- 
sisting only of stronger bonds. Since each cluster is su- 
perfluid by the choice of parameter, the charging energy 
of a cluster of size I is given by A{1) ~ Aq/?. Because 



clusters are independent from each other, the distribu- 
tion of their sizes is simply P{1) = exp{—l/ L), the 
average charging energy on the coarse grained lattice is 
A ~ AqItlL/L. Putting everything together we find: 

Charging energy is thus relevant for a < 0, and irrelevant 
for a > 0, in agreement with the RG flows in Fig. 2. 

The flow towards an insulator for a < implies that 
the superfluid stiffness is suppressed exponentially with 
the system size. But when a = 0+, p{L) ~ 1/lnL where 
L is the system size. This is because the smallest Ji on a 
chain is of order l/L and serves as a lower cutoff for the 
average (1/Ji). Thus in finite but large chains the stiff- 
ness appears to jump at the point a = 0, in analogy with 
KT transition in a pure system. Disorder suppresses this 
jump logarithmically with increasing size. Note that in 
the superfluid regime, a > 0, the collective modes can 
be described by a harmonic chain with off diagonal dis- 
order (see e.g. Ref. [17]). In particular, all finite energy 
phonons in the disordered superfluid are localized [18]. 

In this letter we studied a disordered 1-d 0(2) quantum 
rotor model. Using RSRG analysis, we found a strong- 
randomness fixed point that controls a transition between 
an incompressible Mott glass and a superfluid phase. The 
new phase transition is consistent with KT universality 
class and has a dynamical exponent z = 1. Surprisingly, 
this fixed point lies on the classical axis corresponding 
to vanishing charging energy and uniform distribution 
of Josephson couplings. Unlike other fixed points found 
using RSRG, this one does not belong to the infinite ran- 
domness universality class. Nevertheless, we manage to 
find the physical properties of the model in the critical 
region and in the insulating and superfluid phases. One 
can also use the (universal) solutions of the RSRG flow 
we found to calculate additional thermodnamic quanti- 
ties, as well as finite temperature properties of the sys- 
tem. This we leave for future work. 
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